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1.  Introduction 


The  paper  is  the  second  in  the  series  about  the  h-p  version  of  the 
finite  element  method  for  solving  parabolic  partial  differential  equations. 

In  the  first  paper  [1]  we  discussed  the  case  when  in  the  t  direction 
only  one  element  of  degree  q— »»  was  used. 

In  this  paper  we  analyze  the  case  when  in  both  variables,  the  spatial  and 
time,  the  h-p  version  is  used.  We  are  showing  some  essential  differences 
between  the  p  and  h-p  version.  We  will  keep  the  notation  of  [l],  but  it 
is  not  necessary  to  prerequisite  for  the  present  paper. 


2.1.  The  h-p  version  for  the  initial  value  problem  for  an  ordinary 
differential  equation. 

Let  I  =  (0,T),  I  *  tO.Tl,  t  €  I,  X  =  L2( I )  be  the  usual  space  with  the 

norm 

(2.1.1)  II  u||  x  = 

Let 

V  =  <v  €  C°°(I)  |  v(T)  =  0>, 

where  C°°(I)  is  the  space  of  functions  with  all  continuous  derivatives  on  I. 
For  any  X  >  0  and  v  e  ^  we  define 

(2.1.2)  ||v||Y  =  tt-£  +  AvllX* 

A 

•  Civ  4) 

where  we  denoted  v  =  -rr- •  Let  Y.  be  the  completion  of  C  with  respect  to 

ut  A 

the  norm  IMIY  • 

A 

Remark  2.1.1.  In  [1]  we  have  introduced  in  the  norm 


u2dt 


1/2 


(2. 1.3) 


|v|  =  (lj«^*»Avlx)1/2 


sind  have  shown  (Lemma  2.1  of  [1])  that 


(2. 1.4) 


C1l(v||z  £  |iv|(Y  <  C2l|v||z 
XX  X 


with  >0,  C2  <  +oo  independent  of  X  and  v  but  dependent  on  T. 


On  XxY^  we  define  the  bilinear  form 


(2. 1.3) 


B. 


T 

(u,v)  =  u^-^+Xvjdt. 


Further,  let  F  €  be  a  linear  functional  on  Y  .  We  can  define 


Problem  P^.  For  given  F  €  Y^  find  u^  e  X  such 


that 


(2.  1.4) 


^X ( u0 ’  V )  =  F*v)  v  v  e  Y^. 


This  problem  has  been  analyzed  in  [1]  where,  among  others,  the  unique 
solvability  of  it  has  been  proved.  Moreover,  the  solution  of  the  problem  P 

A 

is  a  weak  solution  of  the  initial  value  problem 


(2.  1.5) 


if 


X  +  Xu  =  f 


u(0)  =  aX, 


F(v)  = 


fvdt  +  av(0) . 


Let  now  k  >  0,  an  integer,  A  =  (t^.tg) 

S  (A)  =  {w  |  w  is  polynomial  of  degree  <  k  on  A) 
§)k(A)  =  {w  e  Sk(A)  |  w(t„)  =  0} . 
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We  will  first  consider  the  following  auxiliary  problems: 

[Irl  k 

for  any  \  >  0,  q  11,  k  =  0,1,...,  c  >  0,  find  u>  e  S  (-1,1)  such  that 


(2.1.6) 


where  is  the  Legendre  polynomial  of  degree  k. 

The  following  lemmas  give  us  the  estimates  of  the  solutions  for  the  above 
problems  which  are  important  in  our  further  analysis. 

r  k  1 

Lemma  2.1.1.  The  problems  (2.1.6)  have  the  unique  solutions  u  for  any 
k  =  0,1,...  .  They  satisfy: 


i) 

II-" 

— —  /> 

/  2 

11  X 

-1,1) 

cv 

'  2k+l 

ii) 

a>[kl(l) 

>  0  and 

|Jk] 

(-1) 

t  S  u[kl(l) 

iii) 

if  k 

^  k2  then 

W[k‘ 

](1) 

<  W[k2](l), 

[k] 

Proof.  Let  us  represent  the  solution  w  of  the  problem  (2.1.6)  by 

k 

PjU). 

J-0 

where  Pj  is  the  Legendre  polynomial  of  degree  J  and  use  the  summation 
formula  (see,  e.g. ,  [2]) 

J 

(2.1.8)  Pj(t)  =  ^(2j-4i-l)PJ_2i_1(t). 

i=0 


(2. 1.7) 


u[kJ(t) 


-I' 


,(k] 


Then  we  get  the  system  of  linear  equations  with  the  nonsingular  matrix 
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(2. 1.9) 


X 

(2k-l ) 

X 

0 

(2k-5) 


X 

(2k-3) 

X 

0 


0  0 

0  0 

X  0 

(2k-5)  x 


1 

X 


1_ 

X 


rR[k]l 

pk 

c 

ftIk) 

Pk-1 

0 

„[k] 

Pk-2 

— 

0 

ft[k] 

Pk-3 

0 

ftfk) 

30 

0 

So,  we  can  see  that  the  representation  (2.1.7)  is  unique  and  the  problem 
(2.1.6)  has  a  unique  solution  for  any  k  =  0,1,2 .  Then 


Ik] 


i) 


u  ,  |cPkl 


2k+l* 


,[k] 


Furthermore,  all  the  coefficients  0 .  ,  j  =  0,1 . k,  are  positive, 

which  immediately  implies  (ii).  Indeed,  from  the  property  of  the  Legendre 
polynomials,  we  have 


|o>tkI  (-1 )  | 


=  i]T(-i)V.k]i  <  £ 

j=0  j=0 


ejk]  1  ~  X0jkl  =  wtk](1)- 


By  analyzing  the  system  (2.1.9)  for  k  and  k+1  we  can  observe  that 


R[k+1]  _  [k]  _  1 

0k+l  -  ^k  ~  X* 


and 


„(k+l]  [k]  .  ,  _ 

^k+l-j  Pk-j  for  J  1,Z* 


PQk+l]  >  0. 


.  ,  k, 


The  above  inequalities  imply  that 

k2  kj 


Is i"21  >  Zsjkl1  lf  k2  >  kr 


k=0 


k=0 
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i.e.  , 


u[klJ(l)  <  if  k1  <  k2. 


Lemma  2 .  1.2.  Let  w^kl3  and  w^23  ^  tlie  solutions  of 


u) 


and 


tki  ] 

-^-  +  c.Xu3kl3  =  c.P.  ,  t  e  (-1,1) 
A  1  1  1 


•[k2] 

T^T'^21'  -=2Pk2’ 


respectively.  Then,  for  kj  <  k2  and  c2/ci  =  <r> 

(2.1.10)  w[kll(l)  <  max(l,<r2k2)^k2l(l). 

Proof .  The  values  u|kl3(l)  and  w2kz3(l)  can  be  expressed  in  the  following 


way  ( see  ( 2 . 1.9)): 


u3kl3(l)  =  - 
1  K  \ 


4k2'(i>  •  i 


J^1  f(i,k  ) 
1+  >  - 2 


i=i  (cix) 


2i 


1  + 


f(i,k2) 


1<2 

I 

i=!  <V> 


2i 


with  some  integer  function  f(i,k)  >  1.  Then 

k? 


4M(i>  -  j 


l  + 


fd.kJ 

I 


1-1  (C2X) 


2i 


1  + 


f(i,k2) 


k2 

Z  M2i 

i=l  (<rcl'° 


max(  1 ,  <r2k2)X 


1  + 


f(i,k2) 


k2 

I 

R  (cix) 


2i 


Now,  by  using  Lemma  2.1.1  with  \— >ACj  and  c  =  c^  we  complet**  the  proof.  ■ 
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Let  now 


(2.2. 1) 


-p  version  of  the  finite  element  method  for  the  problem  P 


4N  ;  0  =  •=  ll  <  (2  <  '  <  lN  '  T 


be  a  subdivision  of  (0,T)  into  time  intervals  I  =  (t  , ,t  ),  n  = 

n  n-1  n 


1.2 . N.  of 

the  length  x^  = 

>  t  -  t  ,, 
n  n-1 

and  d  =  x  /2. 
n  n 

Furthermore,  let 

a  *  (q,.q2.... 

•qN),  q„  2  1. 

an  integer, 

,  n  =  1,2 . N, 

and 

(2.2.3) 

S  =  S3"1 

=  (u  e  X  | 

u|T  e  Sqn_1(I  )> 

1 1  n 

n 

(2.2.4) 

< 

II 

Wo 

IQ 

II 

<v  6  Vx  | 

vl.  e  Sqn(I  )>. 

1 1  n 

n 

The  definition  of  implies  the  continuity  of  v  6  V  and  v(T)  =  0. 

We  define  the  h-p  version  for  the  problem  P :  For  given  F  e  Y'  a 

A  A 

mesh  An>  a  =  (qltq2 . q^)  find  u^  e  S  such  that 


(2.2.5) 


B.(u  ,v)  =  F(v) ,  V  v  €  V. 
a  a 


Remark  2.2.1.  The  problem  (2.2.5)  is  actually  the  set  of  N  problems  to  be 
solved  in  the  succession.  Let  us  denote  for  n  =  1,2, ...,N 

S=  (u  :  I  — »F  |  u  is  a  polynomial  of  degree  <  q  -1}, 
nnn  n  n 

V  =(v  :  I  -»R  |  v  isa  polynomial  of  degree  Sq>,  n=l,...,N-l, 

nnn  n  °  n 

V,  =  (vN  :  In— |  vn  is  a  polynomial  of  degree  <  qn>  v^(t^)  =  0> ,  n  =  N, 
and  represent  the  space  in  the  form  of  the  direct  sum 

(2.2.8)  Vn  “  Vnll0Vi21’  n  =  1 , 2,  .  .  .  ,  N  -  1 ,  V^V^11 


where 


n 


=  (v 


111  €  V  |  V[11(t  )  =  0} 
n  nnn 
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v[2]  =  <vl2]  €  V  |  vI2]  =  c  X  >. 
n  n  n  n  n  n 


where  *n  Is  a  (fixed  but  arbitrary)  polynomial  of  degree  £  q^  such  that 

X_(t  1  }  =  0*  =  1' 

n  n- 1  n  n 

Taking  into  account  the  continuity  of  each  function  v  €  V,  vjT  e  V_  we 

If  “ 


see  that 


dim  V  =  dim  S  =  ,»  q 


N 

=  y,  • 

L~,  n 
n=l 


The  result  is  that  the  problem  (2.2.5)  is  replaced  by  the  following  set  of  N 
problems 


(2.2.7) 


B?(u  .v111)  =  F^v111)  V  -[l1  -  ”[l1 


A  n’  n 


n 


V  -  €  V 

n  n 


where  for  n  =  1 , 2, .... N, 


BX(VVn11} 

Ann 


*[  "n('T-^vin]dt 


with 


(2.2.9) 


F^vJ11)  = 


fvj1]dt  +  avj1](0) 


Fn(vf 1J) 
n 


fvflJdt  + 
n 


fvj^jdt  -  B?_1(u  vj2}),  n  =  2,3, 

n-l  A  n-l  n-1 


.  ,N. 


n-l 


and 


vI2l(tn)  =  vIIi1(t J  v  n  =  1.2 . N-l. 

n  n  n+ 1  n 


Let  us  note  that  the  above  sequence  of  problems  is  independent  of  the 

selection  of  y  . 

*n 

Theorem  2.2. 1 .  i)  Let  u  €  S,  v  e  V,  then 
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(2.2. 10) 


11) 

then 

(2.2. 11) 

ill) 

(2.2. 12) 

Proof.  1 
ii) 

(2.2. 13) 

and 

(2.2. 14) 

we  get 

(2.2. 15) 

and 

(2.2. 16) 

For  given 


|Ba(u.v)|  S  I!u||x||v||y 

X 


If  £  q2  £  .  .  .  £  qN  £  1  and  tj/tj  *  *  for  any  I  S  J  ^  N, 


d.(a,N)  -«  inf  sup  |  B.  (u,  v)  |  2:  1  1/2N_1min(l,«r  (2qi+(  1/2)3 ) 

ueS  veV  A  41 

l|u!|x=i  |v|  SI 

Let  v  €  V,  v  *  0,  then 


sup  |B.  (u,  v)  |  >  0. 
ueS 
llu||x=l 


(2.2.10)  follows  from  the  Schwarz  inequality. 
Denoting  by 


u 


n  X,  n 


u2dt 

n 

n 


1/2 


'  II V  „ 

V 


£.*vl 

A  nj 


1/2 


llu||v  = 


/■  N 

I 

tn=l 


II  u, 


n  X,  n 


1/2 


l|v|| 


1/2 


u  €  S,  u|T  =u  eS,  we  will  construct  v  e  V,  such  that 
ln  n  n 


I 

n 


v  e  V  . 
n  n 
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and  using  v  =  v 

as  the  function  from  the  test 

Let 

qn-i 

(2.2. 17) 

V  n  [k] 
“n  =  L  Vn  ' 

k=l 

where 

qn-i 

(2.2. 18) 

un  =  Z  *VPk> 

k=l 

where  P^  is  the  k-th  Legendre  polynomial  on  I  and  €  S^(I  )  is  the 

unique  solution  of 


0) 


(2.2. 19) 


[k] 

B — +  AufkJ  =  P"  t  €  I  ,  k  =  0,1 . 

A  n  k  n 


By  transforming  these  problems  onto  the  reference  element  (-1,1)  we  can  see 
that 

f  1c  1  k 

ul  J  €  S  (-1,1),  k  =  0,1 . q  -  1 

n  ^n 

satisfy  the  following  equations 
~[k] 

U) 

(2.2.20)  -JL—  +  \dlulk]  =  d*/2P.  ,  t  €  (-1,1),  k-0,1 . q  -  1  • 

.  .1/2  n  n  nk  n 

ACl 

n 

1/2 

The  problems  (2.2.20)  are  of  the  type  (2.1.6)  with  different  A^  =  Adn  and 
1/2 

c  -  d  for  different  n  =  1,2, ...,N.  Thus 
n  n 


(2.2.21) 


B>(u n>w n)  = 

Ann 


u 

n^  A  nj 


dt  =  II  u 


n  X,  n 


and 

(2.2.22) 


|W  liY  n  =  UJY  . 
n  Y^, n  n  x, n 
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Now  let 


(2.2.23) 


u  (t  )  ,  . 

w  -  <.  -  ,  n 

n  n 


fqn] 


where  ulMnJ  is  the  solution  of  (2.2.19)  with  k  =  q  .  Obviously, 
n  n 


(2.2.24) 

and 


w  (t  )  =0,  n  =  1,2 . N, 

n  n 


w  (t  ) 


v  L)  \  t  )  J 


n  n 


since 


B?(u  ,wIqnl)  =  0. 
Ann 


Further 


9n  1  In  1 

(2.2.26)  ivv  - 1£  v£k,(t„»  =  i£ 

k=0  k=0 


r9n-l  Il/pf'*1'-* 

I  I 

4c=0  4?=0 


2JEll(Jkl(t  ))2 

£  n  n 


1/2 


<  d'1/22‘1/2q  l«J.qn](t  )||u  ||  , 

n  n  n  n  n  X,  n 


where  we  used  the  Schwarz  inequality  and  Lemma  2.1.1  iii).  Thus 


(2.2.27) 


l|wj|v  „  S  || u 


l«(t)l  In  1 
n  n  ii  i  qn  j  u 

»  »v  _  a  ii«  !iv  „  +  — m - liw„  "v  „ 

YX*  n  Y^.n  |utqnl{t  )(  n  V 

n  n 

)Hu  iv  d1/2 


u  vj  |  w  v.  u  j  u  _  u  / - =- 

.  ..  ,,  n  Mn  n  n  n  X,n  n  /  2 

S  Un  X, n  ,, - / 


v^|a>l4nJ(t  )| 
n  n 


s  l|u||x  fl+ - — 1  S  2qy2||un||Y  . 

n  X,rH  v^T  J  n  n  X,n 
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c\j  x 


Tkl  [iel 

Since  |w„  (t  ,)|  <  wl  J(t  )  for  k  =  0, 1 . q  (Lemma  2. 1.1  ii) 

n  n- 1  n  n  n 

analogously  as  before 


(2.2.28) 


l«  (t  )l  r _  i 

I  w  (t  )  |  S  1  <i>  ( t  1 )  |  •»  r-n,-n - l<qnJ(t  ,)| 

n  n  l  n  nl  n  n- 1 


(V 


S  d‘1/2v£q ni<qni(tn)l!lu  ,,Y 

n  n  n  n  n  x,  n 


Let  us  now  define  the  following  sequence  of  functions  z e  V  . 


(2.2.29) 

and 


ZN  .  0 


ZN-1  €  VN-1 


be  the  solution  of  the  problem 
z 


2t±,w  -  Wl1  -N-l 

X  NM  '  „'<!«-■  >(t  ^Vl'  "-1' 

N-i  ltN-r 


where  ^  is  defined  by  (2.2.19).  Then 

,WN(tN— i )  I 

(2.2.31)  »z„  |Y  „  ,  *  ln-  I  - BP„  lly  N  - 

N  1  YX’N  1  |W^-ll(tN_1)|  qN-l  X'N  1 


If  "N  5  "N-l 

and  dN-l/dN  ” 

<r„  „  ,  So-,  then  from  Lemma  2.1.2 

N,  N“1 

(2.2.32) 

S  ^  ( t )  1  max(  1,  <r2q,*_1 ) 

and 

(2.2.33) 

^N-l^.N-l  * 

.,1/2  1/2  ,,  2q(i-i+l/2»  „  .. 

2  q^maxd.o-  ^Vx.W 

Further, 

we  gel 


We  put 
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zN-2  6  VN-2 


is  the  solution  of  the  problem 


(2.2.34) 


Since 


ZN-2  .  „  WN-l(tN-2)+ZN-l(tN-2)riN-2 

— +XZN-2 - , - V,’  1  €  X! 

"N-2  ^N-25  N  2 


|zR-l(tM-2J 1  <  ^N-l^N-l51  =  ,wN(tN-l)l 


and  from  (2.2.28) 


(2.2.35)  ,wN-l(tN-2)|  *  dN-l  ^N-l^N-l  1  ( tN-l )  1  l|uN-l!lX, N-l 


while 


(2. 2. 36)  '"N^N-l51  5  dN  V^qN|wNqM  ltMJ  1  I,Un”x,  N’ 


we  have  (with  Lemma  2.1.2) 


(2.2.37)  llzN_2llVN_2  * 


dN-l  ^N-l^N-l  ll(tN-l),,,UN-i(lx,N-l 


,WN-2  2  ^N-2* 1 


dN  V^qN,WNqM  (tN),,|uN,lX.NlJl/2  /  2 

- rr - 1 - cL.  7T- - 


N-2  lN-2 


N-2v  2qN_2+l 


S  21/2^»axa.«'2q"-Z“,/2))(IVl«X,N-l*|lVx.N) 


Hence  we  can  define  the  sequence  z .  €  in  the  recursive  way 


(2.2.38) 


ZN  '  ‘  6  V 


i  .  W  (t  , )+z  (t  ,  ) 

n-l  .  *  n  n-l  n  n-l  „n-2  .  _  „  _ 

-j —  +  Az  =  - T7, - 1 - P  •  tel,  n  =  N,N-  l,N-2 . 2, 

A  n-l  )  qn-2  n_1 

n-l  n-l 

ar.a  «_  w,,  f  O  <5  1  Q  ^  Tk.. 


where  u  is  defined  by  (2.2.19).  Then 

n-l 
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(2.2.39) 


i=n 


Finally,  let  v  e  V  such  that 


(2.2.40) 


v  _  =  w  +  z  for  n  =  1,2 . N. 

1 1  n  n 
n 


Obviously  from  our  construction,  v  is  continuous  and  v(T)  =  0.  Moreover, 


B  (u  ,z  )  =  0  for  all  n  =  1,2,.. .,N. 
Ann 


Thus 


I|U1!IX  1  +  I,U2IIX  ?+---+l|uNlfX  N 
(2.2.41)  sup  |B  (u,v)|  i  -  *  X»N 

v€V  A  (||w1+zj|5  ,  +  .  .  .+l|w„+zj|„  Jl/d 


l|v||  SI 
A 


1  1  Y  ,  1 

A 


N  N  YX,N' 


«Ul«X.r»U2*X,2^--*'UN'x.N 


1/2 


l,Ul l*X,  1  +  '  •  •+I|UNIIX,N 


23/2qJ/2(l|Uil|^  1+. . .  -I  Vlx,N+»ax(l^4qi+1)  [  jjl  Vx,  if) 


N  r  N 


,2>l/2 


and  further 


(2.2.42) 


sup  |B  (u,v)|  £ 
V€V 

l|v||  =1 

A 


N  2 
l  HUJx  n 
n=l  n  X’n 


23^!/2(j1"un'x.„*N(N-1)"“‘'’"-4q,tl)J1lu„«X,„) 

rN  2  i1/2 

U'Vx.J 

.  1/2..  , .  2qi  +  l/2, 

4qj  N  maxd.o-  ) 


1/2 


The  inequality  (2.2.42)  immediately  yields  ii). 


iii)  For  any  veV,  v  =v)T  e  V  ,  we  select  u  e  S,  such  that 

n  'I  n 


n 
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(2.2.43) 


U  =  U  -  =  -v  . 

n  'I  n 

n 


Then  from  continuity  of  v  €  V  and  v(T)  =  0, 

T  N  ^ 

(2.2.44)  B(u,v)  =  u£-^+\vjdt  =  ^  (-v^)  (-  X~ +  *vnjdt 

■*0  n=ri  ’’  J 


N  _  *2  N 


n  r  r  ~  c.  >  n  w  e.  n 

•Z  r11'*  Vndt|*Zf  X^-IZ  it(vn 

n=l^I  -*1  J  n=r  I  n=ri 

n  n  n  n 


)dt 


n  n 

N  „  *2  N 


■If  rdt-5Z[Vn(tn>-''2(t„-l>] 

n=lJI  n*l 


n 

N  _  *2 


■t  rdt‘5[vNttN>-v?(to)] 


n=lJ  I 

n 

N  „  *2 


■  if  rdt-|vi(to)  5  0 

n=lJI 


if  v  *  0. 


Theorem  2.2.1  together  with  Theorem  2.2  of  [1]  and  Theorem  6.2.1  of  [3]  (see 
also  Appendix  of  [1])  yields 

Theorem  2.2.2.  There  is  a  unique  u^  satisfying  (2.2.5).  If  uQ  e  X  is  the 

exact  solution  of  the  problem  P .  and  q,  >  q_  >  .  .  .  £  q.  >  1,  t./t.  S  <r  for 

A  1  <c  N  1  J 

1  i  <  J  S  N,  then 

Nmax(  1,  cr2c*1+^ )j  inf  ||u  -w||  j.. 

J  W€S 


(2.2.45) 


II  un  -  u  | 


a"x 


l+4q 


1/2 


2.3.  Comments. 

If  (2.2.11)  is  optimal,  i.e.,  if  there  is  a  sequence  of  (A^.g^.N^)  such 

that 
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dx  (q^Nj)  S  Cq'I^N*1  min{  1 , <r“(2qi *1+(  1/2) } ) 


with  C  independent  of  i,  then  there  is  a  sequence  of  solutions 
such  that 


lu.-u  || 

i  Si  X 


fnffu,  — w|| '  *  Cqi1^Nimax(l,<r2q‘*1+(1/2)). 


weS 


u1  €  X 


(See  Theorem  2.10  of  (4J.  ) 

The  optimality  with  respect  to  the  exponent  of  q  has  been  proved  in 
(1].  Let  us  now  prove  the  optimality  with  respect  to  the  exponent  of  N. 
Consider  the  case  of  a  uniform  mesh 


(2.3. 1) 
with 

(2.3.2) 

Let  qn  =  q  = 

(2.2.3)  v(t) 


0  -  l0  <  ‘i  < 


<  t„  *  T 


and  v  e  V  be  defined  by  (Figure  2.3.1) 


Figure  2.3.1  v(t)  for  N  =  6. 


5 


Now,  let  ,ii0, . . .  ,0„  ,  €  S  denote  the  basis  function  defined  by 

1  2  N- 1 

fl  tel 


(2.3.4) 


Then  we  have 


(2.3.5)  Bx(0, 


^n(t)  =  i 


n 


n  =  1,2, .  . .  ,N. 


0  til 


„■*>  -  [  ”  - 1-*- 
n 


.  N. 


Any  u  e  S  can  be  written  in  the  form 

N 


(2.3.6) 


then 


(2.3.7) 


Obviously 


u  =  /  c  \1) 
/  .  n*n 


with  ||u||, 


1  V  2 

"  N  n* 


n=l 


n=l 


B»lu';i  *  ^[<-i>n“cn(: 

n=l 


2N-2n+l  X 
NX  2N2 


]]' 


N 


(2.3.8)  sup  |  V*  f(-l)n+1c  [• 
"uHx"1  n=l  1 


2N-2n+l 


NX 


*)]■  • 

L  n=l 


211/2 


2,1/2 


Note  that  there  exists  a  constant  >  0  such  that 

rN2  p,l/2 

(2.3.9)  ||v||Y  >  K2  ^  +  1  +  X  J  . 


Combining  (2.3.8)  and  (2.3.9)  we  see  that  if  X  is  sufficiently  large  in 
comparison  with  N  (e.g. ,  X  £  N),  there  Is  a  constant  C  >  0  such  that 
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(2.3. 10) 


inf  sup  |B^(u, v) |  £  CN 

veV  ueS 

|v|y  -l  Hull x<i 


We  have  to  understand  (2.3.10)  as  an  estimate  which  is  uniform  with 
respect  to  X.  We  have  seen  that  for  given  N,  the  estimate  (2.3.10)  holds 
for  X  sufficiently  large  respectively  to  N,  e.g.,  X  >  N. 

On  the  other  hand  we  can  easily  show  that 


(2.3. 11) 


inf  sup  |Bx(u,v)|  >  C(l+X4)  1. 

veV  U€S 

I  V||  Y  =1  lull  X<1 


In  fact,  given  v  €  V,  select  u  =  -  ^.  Then  u  e  S  and  we  have 


(2.3. 12) 


Bx(u,v)  = 


fTr  2  , 

lV~ 


dt  = 


J 


V  *  F<°> 2 


2 

?*• 


Thus,  using  Remark  2.1.1 
(2.3. 13) 


I  u f|  x  <  c ||  v||  Y 


and 

(2.3. 14) 

(2.3. 15) 

Hence 

(2.3. 16) 


,T 


v  dt  S  C 


v  dt  (because  v(T)  =  0), 


My  *  c 
X 


T  r'  2  „  ' 

| (iPXv  dl 

<S  C(l+X4) 

LJo  x  J 

Jl 

J 


•2 

?*• 


B  (u.v)  >  — 2-  ||v||2 
1+X  X 


which  yields  (2.3.11). 

The  assumption  i  S  . . .  2.  seems  to  be  only  implied  by  the 
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JY 

specific  construction  of  the  sequence  (z  }  ,  in  the  proof  of  Theorem  2.2.1. 

n  n=i 

However,  the  counterexample  considered  below  will  show  that  the  assump¬ 
tion  T1/Tj  *  9  ?or  i  *  J  S  N  is  necessary,  i.e. ,  without  this  condition 
d^(q,N)  can  be  arbitrary  small  (for  some  X).  Let  us  consider  only  two 

elements  of  length  1  and  t,  respectively,  and  q^  =  =  1. 

V 

x 
b 

y 

a 

I  l  +  r  r 

Figure  2. 3.  2 


Thus 

(2.3. 17) 

(2.3. 18) 


v(t) 


'(b-a)t  +  a  0  <  t  S  1 
-  — t  +  b( 1+  -)  1  <  t  <  1+t 

T  T 


u(t)  * 


0  <  t  <  1 

1  <  t  <  1+T. 


Then 
(2.3. 19) 
and 


8(u, v) 


x(a-b)  t  yb  .  ^  fx( a+b)  ^ ybTl 

~  X  *[  2  2  J 


(2.3.20) 

|v||| 

*X 


a2-ab+b2  b2  A  ^2fa2-ab+b2  .  b2r| 

^  [  3  “J' 


Obviously, 


18 


(2.3.21) 

R  =  A  sup  Ix^^  +  yi/r 

tx.y]  2  2 

2  2  . 
x  +y  T=1 

=  A[la2  +  i2ab  +  b2  • 

Let 

(2.3.22) 

inf  A2  j^ia2  +  ^2ab  +  b2 

where  inf 

is  taken  over  a,b  such  that 

(2.3.23) 

A2^2  +  i2ab  +  b2  ^ 

2  JZ  > 1/2 


Q. 


min’ 


Then  Qmjn  is  the  smallest  eigenvalue  of  the  following  eigenvalue 
problem 


(2.3.24) 


AZ  =  QBZ, 


where  Z  =  (a, b)  , 


(2.3.25) 


'1 

1  ' 

‘1 

1  ‘ 

= 

4 

1 

4 

1+T 

B  = 

3 

1 

6 

1+T 

4 

V 

4 

J 

6 

3 

det(A-QB)  =Q2  ^-0^111^  =  0 


and  Qmln— *0  as  r— *0.  Assuming  that  (2.3.23)  holds,  it  is  obvious  that 
there  exists  a  function  0  <  H(t)  <  m,  0  <  t  <  «,  such  that 

(2.3.26)  ||v||2  S  H(t)L+1 

A 

and 

(2.3.27)  sup  |B  (u,v)|  S  K(t)^+R. 

llu||xSl  A  AJ 


Hence,  using  Remark  2.1.1,  we  get 
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d.(x)  =  inf  sup  |B  (u,v)|  5 

|v|Y  -1  Hull X^1  A  ^  X 


Now,  if  X  >  Q  */2H1/2(x)  and  X  >  then 

min 


1/2 


(2.3.28) 


d  (x)  <  2Q1/Z 
X  min 


and  d^  ( ^  j  ( x )  — >  0  as  x— »  0 . 

We  have  once  more  analyzed  here  the  estimate  which  is  uniform  with 
respect  to  X.  Analogously  as  before,  we  could  derive  an  alternative,  not 
uniform  estimate. 

Let  us  now  analyze  in  detail  the  finite  element  method  (2.2.5)  for  the 
problem  P ^  with  q  *»  1,  x  «*  x  *  T/N,  n  =  1,2 . N.  We  get  now 


aX  ♦ 


f(t)(x-t)dt 


(2.3.29) 


W^n  .  X2, 

—  •i-'VVi1  ’  A 


n+1 


f ( t )g  ( t )dt, 


'n-1 


where 


gn(U  - 


— (t-t  )  for  t  ,  <  t  <  t 
x  n  n-1  n 

i(t  ,-t)  for  t  <  t  <  t  . 
x  n+1  n  n+1 


and 

y  =  u( (n-|)x)  =  u  |  _  i  o  m 

n  3  2  3  In(tn)>  H'1’2 . N' 

We  see  that  the  finite  element  method  is  the  well  known  Crank-Nicholson 

difference  scheme.  By  Theorem  2.2.2  we  have 


(2.3.30) 


|uo'Vx  s  CN  inflluo"wi,X’ 

weS 
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where  is  the  exact  solution  of  the  problem  P^.  Obviously,  for  smooth 

2 

function  as  e.g. ,  u^  =  t  we  have 

inf  ||u  -  w|| y  >  CN"1 
weS  U  * 

and  (2.3.30)  suggest  that  the  method  does  not  converge. 

On  the  other  hand,  by  the  classical  finite  difference  approach  we  have 

(2.3.31)  |u0(tn}  "yn'  S 

Letting  PuQ  €  S  be  such  that  Pu^lj  =  uQ(tn)  we  8et 

n 

(2.3.32)  Iu0-uaIx  S  >u0  '  Pu0,X*  ,Pu0  "  ug*X  5  (I)’ 

00 


0N_1  ||Uq  h l  (i) 
00 

CN-2iu<3,lL  (I) 


Hence  we  obtained  the  convergence  in  contrast  to  the  estimate  (2.3.30). 

Let  us  first  note  that  using  (2.3.11)  we  get  instead  of  (2.3.30),  where 
C  is  independent  of  A,  the  estimate 

(2.3.33)  l,u0~VX  1  C(l+A4)inf||u  -w|| 

a  weS 


Further,  we  remark  that  the  estimate  (2.3.31)  assumes  high  smoothness, 
for  example 


Let, 


-A2t 


Then 


uQ(t)  =  aAe 


00 


and  for  A  =  N  we  get  from  (2.3.32) 


,U0_Uglx  S  CN_1a5  =  Cn4> 


while  by  (2.3.30)  we  get  a  bounded  error. 


3.  The  h-p  version  for  the  parabolic  problem. 


3.1.  Preliminaries  and  problem  formulation. 


Let  ft  c  IR  be  a  bounded,  Llpschltz  domain  with  a  piecewise  analytic 
boundary  T.  Let  D  =*  Ixft,  I  =  (0,T).  Then  we  will  consider  the  problem 

3u 


(3. 1.1) 


—  -  Au  =  f  on  D 
o  t 


u  =  0  on  Ixf 


u(0,x)  =  g(x)  on  ft. 


Let,  as  In  [1],  X  =  X(D)  =  L^I.ft^ft))  with 


(3.  1.2) 


fT 

'""x  -  J  '“V 


dt 


(ft) 


and  Y  =  Y(D)  be  the  completion  of 


=  {v  e  C°°(I)  |  v(t,x)  has  for  any  tel  compact  support 
In  ft  and  v(T,x)  =  0} 


In  the  norm 


(3.  1.3) 


.  •  3v 

where  v  = 

o  t 


,T 


civi"'  +  \m.  )dt 

,  H  (ft)  ft  (ft) 


On  XxY  we  consider  the  bilinear  form 


(3.1.4) 


B(u, v)  = 


(~uv+7u7v)dxdt 


CT 


ft 


and  the  problem  P:  Given  F  e  Y' ,  find  u^  e  X  such  that 
(3. 1.5)  B(u, v)  =  F(v)  V  v  e  Y. 

In  (1)  It  has  been  shown  that  problem  P  has  a  unique  solution  for  any 


22 


F  €  Y'  and  the  solution  u^  of  the  problem  P  Is  a  weak  solution  of  (3.1.1) 
with 


F(v) 


g(x)v(0, x)dx  + 

n 


|  f (t,x)v(t,x)dxdt. 

■Vo 


3.2.  The  semidiscrete  problem.  Discretization  In  x. 

ol 

Let  R  c  H  (£1)  be  a  finite  dimensional  subspace  of  functions  and 
S  =  (u  €  X  |  u(t,x)  e  R  V  t  6  I>, 

V  =  (v  e  Y  |  v(t,x)  €  R  V  t  e  I>. 

The  semidiscrete  approximation  of  the  solution  for  the  problem  P  is  defined  as 
the  function  ug  e  S  satisfying 

(3.2.1)  B(u  , v )  =  F(v)  V  v  e  V. 

s 

The  following  estimate  of  the  error  for  this  semidiscretization  has  been 
derived  in  [1], 

(3.2.2)  l|uO~usl!X  "  c(*(R))"1irrfllVusllX' 

W€S 

under  the  assumption  that  the  space  R  has  the  property  X,  l.e.,  there  is 
a  number  X(R)  <  1  such  that 

(3.2.3)  ||u||  >  K(R)||u|| 

(n)  H  1  (£1) 

holds  for  any  u  e  R  with 
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Moreover,  it  has  been  proven  in  [1}  without  assuming  the  property 


(3.2.4) 


■Wx  *  C|'VPOuO!lX’ 


where  Pq  denotes  the  Lg-orthogonal  projection  of  X  onto  R. 
Assuming  that  R  is  (£, tj, u)-regular,  i.e., 


II  u| 


£(R)  =  sup 


ft*(n) 


ueR  l|u,1L2(n) 


<  +00 


T)(R)  =  sup  ||U  -  P1uilL  (n)  <  +00 

ueft^fl)  2 

Hull  i  SI 

ft*(n) 


u(R)  =  sup  ||u-P0u||L  (Q)  <  +oo, 

ueft*(n)  2 

l|ull  1  SI 
ft  (fl) 


where  PQ,  respectively  P^,  denotes  the  L^-,  respectively,  HJ 
operator  of  ft*  (II)  onto  R,  we  have  proved 


(3.2.5) 


H(R)  >  (1+?(R)t}(R))-1. 


Further  we  have  shown  that  for  any  u  e  ft* (ft),  ||u||  ,  =  1 

ft*(Q) 


(3.2.6) 


|u-Pnu||  ,  <  (1+t)(R)+v(R))?(R). 

0  ft*(n) 


3.3.  Numerical  examples. 

Let  us  consider  the  following  problem: 

u  -  u"  =  0  for  (t,x)  e  D  =  (0,l)x(0,l) 
(3.3.1)  u'(t,0)  =  0,  u( t , 1 )  =0  for  t  e  (0,1) 


X,  that 


-projection 
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u(x,0)  =  g(x)  for  x  €  (0,1), 


where  u  =  u'  =  u"  =  — Due  to  the  symmetry,  the  above  problem  is 

equivalent  to  the  problem  on  D  =  (0, 1 ) x ( — 1 , 1),  homogeneous  Dirichlet 

boundary  conditions,  and  a  symmetric  solution. 

2 

We  will  address  the  cases  g^(x)  =  1-x  and  (x)  =  1  -  x.  In  the 

first  case,  the  solution  has  singularity  in  the  point  x  =  1,  t  =  0,  while  in 
the  second  case  the  singularity  is  located  in  the  point  x  =  0,  t  =  0. 

First,  we  consider  the  case  of  one  single  element  in  time  and  three  dif¬ 
ferent  meshes  in  space  with  discretization  in  space  only.  The  first  mesh  (a) 
consists  of  one  element  only.  The  mesh  (b)  and  (c)  is  composed  by  two  ele¬ 
ments  with  the  nodal  point  in  x  =  0.05  and  x  =  0.95,  respectively,  and  the 
space  S  polynomials  of  degree  p  on  each  element.  Obviously  for  g 


>1 


inf  tl u  —  w|| 
weS 


is  essentially  the  same  for  the  meshes  (a)  and  (b)  and  the  asymptotic  rate  is 
the  same  because  the  refinement  is  at  the  place  where  there  is  no  singularity. 
In  the  case  of  mesh  (c)  we  expect  essentially  two  phases,  the  first  one  when 
the  rate  is  exponential  (in  p)  and  the  second  one  for  sufficient  high  p 
when  the  rate  is  algebraic,  the  saune  as  for  meshes  (a)  amd  (b)  (for  more  see 
(51). 

On  the  other  hand  tt(R)  is  different  for  the  meshes  (a)  atnd  (b)  and 
H(R)  is  the  same  for  the  meshes  (b)  amd  (c).  In  fact,  (3.2.5)  yields  that 


m R)  > 


h 

max 

h  , 
min 


where  h  ,  respectively  h  ,  is  the  length  of  the  maximal,  respectively 
max  min 

minimal,  element.  Figure  3.3.1  shows  that  the  effect  of  H(R)  does  not 
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appear  In  the  computations. 

For  the  mesh  (a)  we  have  two  estimates  (3.2.2)  and  (3.2.4).  The  estimate 

(3.2.4)  leads  to  the  rate  0(p  ^ e)  as  shown  in  [1],  while  the  estimate 

(3.2.2)  gives  the  rate  0(p  ),  e  >  0,  arbitrary.  The  fact  that  the 

meshes  (a)  and  (b)  give  practically  identical  results  shows  that  the  factor 

H(R)  does  not  influence  the  results.  Further,  the  case  of  the  mesh  (c)  shows 

exactly  what  had  to  be  expected,  namely  a  concave  curve  (which  would  straight- 

5 

en  for  higher  p).  The  theoretical  slope  p  based  on  3.37a  [1]  is  shown 

in  Figure  3.3.1,  too. 
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SPACE  DEGREE  p 
3  4  6  8  10  12  14  16 


Figure  3.3.1.  The  relative  error  of  the  semidiscrete  method  (discretizat ion ) 
in  x)  vs.  the  space  degree  p  in  the  log  log  scale  for  the 
problem  (3.3.1)  with  gj(x)  =  1  -  x2  (slope  based  on  3.37  a 
[1]). 

(a)  one  space  element  h  =  1. 

(b)  two  space  elements  ht  =  0.05,  h2  =  0.95 

(c)  two  space  elements  hi  =  0.95,  hg  *  0.05. 
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RELATIVE  ERROR 


Figure  3.3.2  shows  the  analogous  results  for  gg(x)  *  1  "  x-  Now  the 

singularity  is  in  x  *  0,  t  *  0  and  hence  the  mesh  (b)  is  better  than  mesh 

“2  S 

(c)  in  contrast  with  the  previous  case.  The  theoretical  slope  p  is 

displayed  in  Figure  3.3.2,  also. 

Hence  we  can  conclude  that  the  performance  of  the  method  Is  not  sensitive 
to  the  mesh  (at  least  in  the  one-dimensional  case). 


SPACE  DEGREE  p 

2  3  4  6  8  10  12  14  16 


"1 - 1 - H 

g2 (x)  =  I  - x 


Figure  3.3.2.  The  relative  error  of  the  semidiscrete  method  (discretization 
in  x)  vs.  the  space  degree  p  in  the  log  log  scale  for  the 
problem  (3.3.1)  with  g2(x)  =  1-x  (slope  based  on  section 
3.4  (11). 

(a)  one  space  element  h  =  1. 

(b)  two  space  elements  hj  =  0.05,  h2  =  0.95 

(c)  two  space  elements  hj  =  0.95,  h2  *  0.05. 
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3.4.  The  semi discrete  problem.  Discretization  In  t. 


Let  now  (as  in  Section  2.2) 


(3.4. 1) 


4N  :  0  ■  *0  <  ll  <  C2  <  •  •  <  CN  ■  T 


be  a  subdivision  of  (0,T)  onto  time  intervals  I  =  (t  , ,t  )  of  the  length 

n  n- 1  n 

x  *  t  - 1  Furthermore,  let 

n  n  n-1 

(3.4.2)  g  =  C q  1  * q2 . qN^’  qn  ~  1’  ^  inte8er*  n  =  1.2 . N 

and 

(3.4.3)  T  =  (u  e  L2(I)  |  u I _  e  Sqn~ 1 ( I  )> 

a-l  1 1  n 

*  n 


(3.4.4) 

Define  now 

(3.4.5) 

(3.4.6) 


f*  =  (v  e  C°U)  I  v|T  €  Sqn(I  )  and  v(T)  «  0} . 
g  'In 

n 


s  =  t  xft^n) 
a  a-i 

v  »  f’xft^n). 
a  a 


For  any  u  e  S^,  respectively  v  e  V  ,  we  have 


(3.4.7) 


u( 


t,x)  =  ^ai(t)u1  (x). 


i  =  l 


where  u^  are  the  eigenfunctions  of  the  corresponding  elliptic  eigenvalue 
problem  introduced  in  [1],  with 


a 


i 


e  T 


a- 


1 


and 


(3.4.8) 

with 


v(t,x)  =  ^(thyx), 
i=l 
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Defining  the  bilinear  form 

(3.4.9)  B(u, v) 


e  *a 


B(u, v)  on  S  xV  by 
9  9 


-f[i> 

J0*-1=1 


C-aiPi +  XiaiPi) 


dt. 


where  A^  is  the  eigenvalue  corresponding  to  u^,  i  =  1,2 .  we  can  use 

Theorem  2.2.1  to  prove 

Theorem  3. 4. 1.  Let  u^  be  the  solution  of  the  problem  P  and  us  its  semi- 
discrete  solution  (discretization  in  t).  Then,  if  q  £  q  £  ...  £  q  £  1 
and  Ti/Tj  ^  *  f°r  1  M  <  J  i  N  then 

||  u  -  u  ||  <  CqJ/2N  max(l,<r2qi+(1/2))inf  ||u  -w|| 

°  weS 

9 


3.5.  Numerical  examples. 

Let  us  consider  first  the  case  when  the  initial  function  g^(x)  =  cosg~ 
and  the  semidiscrete  (discretization  in  time)  is  used.  In  this  case  the 
problem  essentially  reduces  to  the  ordinary  differential  equation  because 
g^(x)  is  an  eigenfunction.  Figure  3.5.1  presents  the  results  for  different 
meshes.  Theorem  3.4.1  gives  the  error  estimate  which  strongly  depends  on  <r, 
Figure  3.5.1  shows  that  the  factor  depending  on  c r  is  not  essential,  mainly 
because  the  eigenvalue  A  is  not  too  large  and  the  major  effect,  as  before, 
arises  from  the  approximation.  We  mention  that  we  could  give  for  this  case 
another  estimate  where  the  right  hand  side  depends  on  A. 
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TIME  DEGREE  q 

I  2  3  4  5  6  8 


Figure  3.5.1.  The  relative  error  of  the  semidiscrete  method  (discretization 
in  t)  vs.  the  time  degree  q  in  the  log  log  scale  for  the 
problem  (3.1.1)  with  g3(x)  =  cos^  and  the  following  meshes: 

(a)  Ti  =  0.1,  t2  =  0.2,  t3  =  0.3,  r4  *  0.4,  <r  <  1 

(b)  xt  =  0.1,  x2  =  0.4,  T3  =  0.3,  x4  =  0.2,  cr  =  2 

(c)  rx  =  10/71,  t2  =  50/71,  r3  =  1/71,  x4  =  10/71,  <r  *  50. 

In  the  other  example  (Figure  3.5.2)  we  did  not  observe  the  dependence  on  <r 
(although  for  some  solutions  such  dependence  will  occur,  see  (2.3.28)). 
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TIME  DEGREE  q 

I  2  3  4  6  8  10  12  16 


Figure  3.5.2.  The  relative  error  of  the  semidiscrete  method  (discretization 
in  t)  vs.  the  time  degree  q  in  the  log  log  scale  for  the 
problem  (3.1.1)  with  gi(x)  »  1  -  x2  and  the  following  meshes: 

(a)  n  »  1/85,  x2  *  4/85,  t3  *  16/85,  x4  *  64/85,  <r  <  1 

(b)  Ti  =  1/85,  x2  =  64/85,  x3  =  16/85,  x*  =  4/85,  or  =  16 

(c)  Xi  =  4/85,  x2  =  64/85,  x3  =  1/85,  x4  =  16/85,  «r  =  64. 


3.6.  The  complete  discretization. 

Let  us  define 

(3.6.1)  S  =  S(R,a,A)  =  T^jXR 

(3.6.2)  V  =  V(R,q,A)  =  f^xR. 

Combining  the  results  of  previous  sections  we  get 

Theorem  3.6.1.  Let  un  €  X  be  the  solution  of  the  problem  P  and  u  €  S  be 
the  finite  element  solution 
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(3.6.3)  B(u  , v )  =  F(v)  V  v  e  V. 

s 

Then,  if  the  space  R  has  the  property  X,  q.  >  q_  >  .  .  .  >  q.  >  1,  and 

Id  N 

x^/Tj  5  <r  for  1  <  i  <  j  £  N,  then 

(3.6.4)  #us-u0llx  *  CqJ/2N(mR))"1max(l,o-2qi+U/2))inf||u0- w|ls 

weS 

and 

Theorem  3.6.3.  Let  e  X  be  the  solution  of  the  problem  P  and  ug  e  S  be 
the  solution  of  (3.6.3).  If  denotes  the  error  of  the  semidiscrete  method 

in  (3.4.9)  (discretization  in  t)  and  is  the  error  of  the  semidiscrete 

method  (3.2.1)  (discretization  in  x)  and  R  is  (£, tj, lO-regular,  then 

(3.6.5)  ||us  -  UqII x  *  C(ct(  1  +  r?(R)+p(R))?(R)+e2). 


3.7.  Numerical  examples. 

The  h-p  version  of  the  finite  element  method  gives  large  freedom  to 
select  elements  in  the  space  and  time  variables.  This  flexibility  can  be 
employed  in  various  ways,  for  example,  in  connection  with  adaptive  approaches. 
We  also  have  seen  that  for  q  =  1,  the  method  coincides  with  the  Crank- 
Nicholson  method  and  hence  the  h-pversion  in  t  cam  be  implemented  as  a  spe¬ 
cial  solver  for  the  stiff  ODE’s  arising  From  the  solution  of  the  semidiscrete 
method  discussed  in  Section  3.2.  Various  aspects  of  this  feature  will  be 
addressed  in  the  forthcoming  papers  in  this  series.  Here  we  will  present  some 
illustrative  results  related  to  the  question  of  the  optimal  relation  between 
p,q  and  a  mesh  in  one  dimensional  setting  we  used  earlier. 

We  will  assume  that  in  space  (i.e.,  x-variable)  only  one  element  of 
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degree  p  is  used,  while  in  the  time  variable  we  use  N  elements  of  degree 
q.  The  shape  functions  are  integrals  of  Legendre  polynomials  (in  x),  and  in 
t  we  use  Legendre  polynomials  as  trial  shape  functions  and  integrals  of 
Legendre  polynomials  and  linear  ones  as  test  functions.  Using  the  band  solver 
we  need 

(3.7.1)  W  =  N*  W 
arithmetic  operations,  where 

3 

p  for  2q  i  p 

(3.7.2)  W  ~  • 

_4q  p  for  2q  <  p. 

We  consider  the  problem  (3.3.1)  for  g(x)  =  gg(x)  =  1-x  and  g(x)  = 
ftx 

g^(x)  =  cosg—  as  ttie  representatives  of  the  solutions  for  unsmooth  and  smooth 
initial  data  problems. 

In  the  case  of  g(x)  =  g^(x)  we  have  used  the  radical  mesh  in  t,  i.e. 

(3.7.31  tn-  (a]\  n-  0,1,2 . N.  with  7  =  3. 

and  for  g(x)  =  g^(x)  we  have  used  the  uniform  mesh 

(3.7.4)  t  =  £,  n  =  0,  1,2 . N. 

n  N 

In  the  following  figures  we  present  the  accuracy  in  dependence  on  (p.q) 
in  the  scale  logtrelative  error)  vs.  a  with  qa  =  p.  We  also  show  the  work 
needed  for  computation  based  on  (3.7.1  -  3.7.2). 

For  low  q,  the  error  is  governed  by  the  time  integration,  i.e.,  we  deal 
here  essentially  with  the  case  of  semidiscrete  method  with  time  discretization 
(as  discussed  in  section  3.4). 

In  contrast,  high  q’s  essentially  show  the  performance  of  the  semidis¬ 
crete  method  with  discretization  in  space  only  (as  discussed  in  section  3.2). 
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In  Figures  3.7.1  -  3.7.3  we  present  the  results  for  g(x)  =  g^Cx)  ~ 

1-x.  We  see  that  for  any  required  accuracy  the  minimal  computational  work 
is  needed  when  q  is  selected  low  and  N  large. 

Note  that  the  numbers  in  parentheses  indicate  the  value  of  q  used  and 

4 

e.g.  the  value  5.3E4  denotes  that  the  approximate  work  W  «  5.3*10  has 
been  needed  to  compute  the  finite  element  solution. 

In  Figures  3.7.4  -  3.7.6,  we  show  the  analogous  results  for  g(x)  = 
g^(x)  =  cos^.  Here  we  see,  in  contrast  to  the  previous  case,  that  the  best 
is  to  select  one  time  element  and  large  q. 

In  both  cases  we  used  the  same  p  and  q  in  all  time  intervals.  The 
results  show  that  the  flexibility  of  the  method,  when  properly  employed,  leads 
to  the  large  increase  of  computational  effect ivity.  Various  adaptive 
approaches  here  will  be  very  effective  tools  for  such  an  optimal  choice. 

These  aspects  will  be  addressed  in  the  forthcoming  paper. 


i 
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dOdd3  3AllV“!3d 


Figure  3.7.2.  The  performance  of  the  h-p  version  (N  =  6)  for  the 

relationship  p  *  qa.  The  initial  function  is  g2(x)  =  1  -  x. 


36 


Q2^x^  =  J“x  Radical  mesh 


dOUU3  3AI1V13& 


Figure  3.7.3.  The  performance  of  the  h-p  version  (N  =  12)  for  the 

relationship  p  =qa.  The  initial  function  is  g  (x)  =  1-x 

b 
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cm  ^  cd  co  o  cm 

*  O  'o  'o  o  'o  *C3 


d0dd3  3AllV13d 

Figure  3.7.4.  The  performance  of  the  p-version  (N  =  1)  for  the  relation¬ 
ship  p  =  qa.  The  initial  function  is  g^(x)  *  cos™. 
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